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ABSTRACT 


Teeth provide meaningful clues of an individual. The growth of the teeth is 
correlated with the individual age. This correlation is widely used to estimate 
age of an individual in applications like conducting forensic odontology, 
immigration, and differentiating juveniles and adolescents. Current forensic 
dentistry largely depends on laborious investigation process that is performed 
manually and can be influenced by human factors like fatigue and 
inconsistency. Digital panoramic radiograph dental images allow 
noninvasive and automatic investigation to be performed. This paper presents 
analyses on third molar tooth segmentation for the population in Malaysia, 
ranging from persons age of 5 years old to 23 years old. Two segmentation 
techniques: gradient inverse coefficient of variation with dynamic 
programming (DP-GICOV) and Chan-Vese (CV) were employed and 
compared. Results demonstrated that the accuracy of DP-GICOV and CV 


Third molar tooth were 95.3%, and 81.6%, respectively. 
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1, INTRODUCTION 

The development of a human body can be measured with the biological age as its indicator. The 
biological age of an individual may not be the same with his or her chronological age, due to development 
variations of the body. Development of dental, bone, menarche, voice change, height, and weight can be used 
as variables to measure the biological age of the human body. The development of dental is more reliable 
than other variables to indicate the biological maturity in children and also more resistant than bones in the 
taphonomic process [1]. The measurement of the biological age based on dental can be associated with 
chronological age of an individual. The effort to associate the dental age with chronological age can be traced 
back to the year 1837 when an English dentist, Edwin Saunders presented a pamphlet entitled “The Teeth a 
Test of Age” in the British Parliament. His writing was related to child labor in factories in the United 
Kingdom [2]. Dental age estimation is an important information for forensic odontologist in order to assist 
government agencies and authorities to solve several cases, such as, determining the legal age for child in 
employment and criminal law or their eligibility to receive social benefits, processing undocumented 
immigrants, narrowing the search and identifying unknown victims. 

The third molar is the most suitable teeth for dental age estimation in forensic investigation [3]. The 
third molar is related with the chronological age because it is the only tooth that still growing until late 
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adolescence, which can be regarded as legal age in the court of law. For an example, the legal age in the 
United States (US) is 18 years of age. Other countries may have different legal ages. Even so, the third molar 
growth display varying rates across races and sexes. In the US, the development of the third molar in 
American Hispanic is about one-half year earlier than American White, while for the American Black, their 
third molar develop one-half year earlier than American Hispanic [3]. Comparison between nine countries: 
Belgium, China, Japan, Korea, Poland, Thailand, Turkey, Saudi Arabia, and India also demonstrated 
significant differences in the development of the third molar. In addition, their finding also revealed that at all 
ages, the third molar growth of males is ahead of females [4]. Similar findings are also reported in China [5], 
Brazil [6], and Croatia [7]. All these studies suggest that the development of the third molar is influenced by 
ethnicity, hereditary, genders, nutrition, disease, medical treatments, habits, addictions, occupation, and other 
local effects. However, the most important factor is ethnicity and genders [3, 7]. 

There are methods to estimate the age from dental records [1], such as, morphological, biochemical, 
and radiological methods. Morphological method is the teeth assessment whereby the teeth need to be 
extracted from the human body. The morphological method examine the relationship between the 
chronological age and the changes in the dental tissue. These changes are attrition or tooth wear, periodontitis 
or gum diseases, secondary dentin, cementum apposition, root resorption, and transparency of the root [8, 9]. 
Biochemical method involves amino acid racemization, an organic chemistry process that occur in dental 
biopsy to estimate the age of an individual [10, 11]. Biochemical is a multidisciplinary approach involving 
forensic odontology, biochemistry, histology and metric analysis [12]. Unlike the morphological method, 
biochemical method can be performed without extracting teeth. Radiology method uses panoramic x-ray 
dental images, which allow teeth investigation without removing the teeth form the body [5]. 

Radiology is an effective method because it is a noninvasive procedure and need not to remove the 
teeth from the body. The panoramic radiography images that are stored digitally can be analyzed using 
advanced image processing and other numerical techniques. Deep convolution neural network technique is 
used to perform teeth classification [13-16]. This automatic classification helps inexperience dentist to 
overcome the laborious task of manually recording the dental chart for corpses. Augmented reality and image 
registration techniques are used in dental surgery [17]. Image subtraction based segmentation is used to 
quantify leakage in dental restorations [18]. Multiple morphological gradient technique that consist of image 
dilation and erosion is employed to enhance the identification of proximal caries or tooth decay [19]. Volume 
rendering and texture analysis techniques are used to detect inferior alveolar nerve canal in dental implant 
surgery [20]. However, most of the image processing techniques are only focusing on teeth classification and 
did not relate the growth of the teeth with the chronological age of persons or victims. 

In this paper, the segmentation of the third molar teeth is discussed. DP-GICOV [21] is employed in 
the segmentation and its results are compared with other segmentation technique, CV algorithm [22]. The 
segmentation is performed for a Malaysian dataset ranging on persons aged 5 years-old to 23 years-old. 
Results of segmentation are discussed with respect to the age ranges. 


2. RESEARCH METHOD 
2.1. Panoramic radiography dental image 

Digital images allow noninvasive analysis to be performed. In the forensic odontology, an example 
of a digital panoramic radiography dental image of a person is shown in Figure |. Out of the many teeth, 
third molar or wisdom tooth is the most reliable indicator for age estimation. In a normal person, there are 
four third molar teeth as shown in four boxes in Figure 1. Two third molar teeth that are located at the upper 
dental jaw are not considered in this study because they tend to overlapped with their neighbors (i.e. second 
molar) as the person grow older. This situation may cause a complication in the image segmentation process. 
Only the third molar teeth at the lower dental jaw are considered because they are visibly apparent 
throughout the chronological development of a person and not affected by the overlapping problem by other 
teeth. Here, either one the teeth at the left or mght can be used in the image processing of 
forensic odontology. 
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Figure 1. A panoramic radiography dental image 


Four highlighted boxes show four wisdom teeth or third molar teeth. Only one of the teeth at the 
lower dental jaw is used in the image processing of forensic odontology. 


2.2. Dynamic programming-gradient inverse coefficient of variation 

Gradient inverse coefficient of variation (GICOV) measures the ratio of the mean and standard 
deviation of directional image derivatives over an entire closed contour fitted using active contour techniques 
to a tooth. GICOV operates by generating a number of radial lines originating from a centroid of the third 
molar tooth, which can be determined by users. Figure 2 shows a number of radial lines with different 
orientation from 0° to 360° angles. On each radial line, a number of sampling point is generated. Here, there 
are M radial lines, and N sampling points on each radial line. 


(c,+ R cosG c, + R sin@) 





contour 





Figure 2. Generating M radial lines and N Figure 3. Converting Cartesian coordinate to 
sampling point on each radial line. The polar coordinate 
centroid is manually defined by users 


In order to focus the attention to the sampling points on every radial lines, the Cartesian coordinate 
of the contour points are converted to the polar coordinate as (c, + R cos @,c, + R sin @) and shown in 
Figure 3. The gradient of the derivative images is computed as in (1). 
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g= fl +I (1) 


y 


where J, = O1/Ox and J y= /Oy. I is the image with polar coordinates. Further, directional information 


of the gradient is computed in (2) 
g =§%-sign(1, cosO+ I, sin 0) (2) 


where g is the angle of the radial line taken from a centroid of a two-dimensional closed contour. This 
centroid can be used as a point of initialization for the closed contouring method. The directional derivatives 
are the outward normal directions on contour points [21 ]. 

Figure 4(a) shows intensity profiles taken from point A to point B. This line denotes a radial line 
whereby point A is used as a centroid. Two plots from point A to B are drawn in Figure 4b, which shows the 
intensity plot of original image and the intensity plot of directional image gradient. Along the line from point 
A to B, dynamic programming (DP) [21, 23] is employed to determine a boundary point of the tooth, vj. It is 
observed that the detection of the point v; is easier by using the intensity of directional image gradient rather 
than the intensity of the original image. 
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Figure 4. (a) A radial line generating from centroid, A to point B. vj is the boundary point of the tooth. (b) 
Intensity profiles for original image and directional image gradient taken from the centroid point, 
A to point B 


GICOV score, G is the ratio of the mean and the variance of directional gradient, g(u) at a point uw 
on the first radial line through, g(uy) at a point uy on the N-th radial line as given in (3). 


(mean| g (uy), 8 (Uy )) 


G(uy,,...,4y ) = -—— (3) 


var| g (u, — & (uy )] 


In order to maximize (3), this equation is rewritten as given in (4) 


(haw) (a 
Yorn) -(Yy di alu) +8 


G 


where s is the noise variance. The cost function of the active contour is given in (5). 


EVs) HE, (v,,v,)+E, (v5.3 ) 
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(5) 


Bulletin of Electr Eng and Inf, Vol. 8, No. 1, March 2019 : 253 — 260 


Bulletin of Electr Eng and Inf ISSN: 2302-9285 O 257 


(5) is computed using dynamic programming. Each cost function is defined in (6). 


pny fH) Peale) if D(v,.%1) SB (6) 


CO otherwise 


where function D(x) measures the distance between point v; to point visi, and D is the distance threshold. 
Figure 5 shows the generation of the active contour from a linkage of points v1, v2, ... Vi,... Vl, and 
vu. Here, the combination of DP and GICOV is called DP-GICOV. 





Figure 5. The dashed lines are the radial lines, and each radial line has a point, v; indicating the boundary 
point of the tooth 


All the points, vi, v2, ... Vi... Vwi, and vy are connecting with active contour using dynamic 
programming method. 


2.3. Datasets 

In order to conduct the current study, panoramic radiography dental images of living persons aged 
between 5 to 23 years old with no obvious developmental abnormalities are used. The sample are obtained 
from the Faculty of Dentistry, Universiti Teknologi Mara (UiTM). The third molar teeth either on the left or 
right of the lower dental jaw are analysed. Third molar of the upper dental jaw are not considered. Both 
erupted and non-erupted third molars are included in this study, as long as the third molar itself is 
distinguishable. The panoramic radiography dental images with the following anomalies are filtered out from 
this study: subjects with unknown dates of birth, subjects with missing third molar, subjects whom images 
could not be interpreted properly, and subjects with rotated third molar with large areas of enamel overlapped 
with neighbouring teeth. There are 650 images used in this study. Table | illustrate the distribution age of 
the samples. 


Table 1. Distribution of the samples based on age group 


Group Age No. of Samples 
Group | 5-7 years old 138 
Group 2 8—10 years old 72 
Group 3 11-13 years old 68 
Group 4 14-16 years old 93 
Group 5 17-19 years old 97 
Group 6 20-23 years old 182 

Total: 650 


Figure 6 shows samples of third molar teeth for the different age groups in order to demonstrate the 
morphological development of the third molar tooth. The stages of the tooth morphology starts with the 
appearance of the cusps tips, the merging of the cusps tips to form a well-defined coronal morphology, the 
formation of the crown, the appearance of the pulp chamber in a trapezoidal form, and the formation of the 
inter-radicular bifurcation. At this stage, the length of the root is less than the crown length. Then, the 
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development of the third molar continues with the expansion of the root length twice than the length of the 
crown. At this stage, the roots look like a funnel-shape ending. The root walls appear to be parallel, but tooth 
apices remain open. The development of the third molar concludes when apical ends of the roots are 
completely closed and periodontal ligament has a uniform width around the root. As shown in this figure, the 
progressive morphological growth of the third molar is correlated with the chronological age of an individual 
[3]. Images Figure 6 pose image processing challenges because contrast of a tooth is not uniform. At the 
crown part, the contrast is high between the tooth and the gum, but at the roots part, the contrast is relatively 
low. Figure 6f poses more difficult situation because the contrast of all the part of the tooth and the gum is 
very low. In addition to that, all the images appeared to be grainy, which might affect the accuracy of image 
segmentation. 





(a) (b) (Cc) (d) (e) (f) 
Figure 6. Samples of third molar teeth for persons aged (a) 5, (b) 10, (c) 13, (d) 15, (e) 18, 
and (f) 20 years old 


3. RESULTS AND ANALYSIS 

Figure 7 shows results of DP-GICOV contour using different number of radial lines applied on a 
third molar tooth of a 10 year-old person. DP-GICOV requires suitable number of radial lines. If the number 
of radial lines are too small, the active contour of DP-GICOV is not able to track the boundary of the tooth. 
Using 300 lines allow the DP-GICOV to closely track the boundary. However some parts of the tooth, 
especially the roots of the tooth are not included in the contouring because the intensity contrast between the 
roots of the tooth and the gum is too low. The more radial lines are used, the slower the computational time 
for the DP-GICOV contour. 





(a) (b) (Cc) (d) (e) 
Figure 7. Results of DP-GICOV contour using different number of radial lines. (a) 30 radial lines, (b) 70 
lines, (c) 150, (d) 200, and (e) 300 lines 


Figure 8 shows results of segmentation using DP-GICOV and CV for samples of third molar teeth 
for different age groups. The segmentation accuracy of a sample from Group | (5 years old) using DP- 
GICOV is 94.2% and CV is 87.5%. Group 2 (10 years old)—-DP-GICOV=96.1%, CV=89.5%. Group 3 (13 
years old)—-DP-GICOV=93.8%, CV=78.8%. Group 4 (15 years old)-DP-GICOV=93.5%, CV=80.0%. Group 
5 (18 years old)-DP-GICOV=95.4%, CV=68.0%. Group 6 (20 years old)—DP-GICOV=93.5%, CV=72.8%. 
CV shows poor segmentation results whereby in all the images, this technique is heavily underestimated the 
boundary of the third molar. CV is severely influenced by the low contrast of the roots of the tooth and the 
gum, therefore, this part is not successfully segmented. Moreover, the effect of the grainy image also give 
rise to the inefficiency of the CV technique. Oppositely, DP-GICOV demonstrate close estimation of the 
boundary of the third molar. DP-GICOV is able to segment the roots of the tooth but rather inefficient to 
separate those roots. DP-GICOV demonstrates robustness when dealing with grainy images. 
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(a) (b) (Cc) (d) (e) 
Figure 8. Segmentation results on the third molar tooth on persons aged (a) 5, (b) 10, (c) 13, (d) 15, (e) 18, 
and (f) 20 years old. Red, blue and yellow contours refer to the boundaries of ground truth, CV and DP- 
GICOV, respectively 


Figure 9 shows the histograms of segmentation accuracy of DP-GICOV and CV for 650 third molar 
teeth samples that are divided into six age groups ranging from persons aged 5 years old to 23 years old. The 
morphology structure of the teeth changes as they undergo the teeth grow that is according to the 
chronological age of human being. From the histograms, there are two patterns between the accuracy of DP- 
GICOV and CV. DP-GICOV demonstrates a consistent pattern for all age groups. DP-GICOV is able to 
perform segmentation for approximately full shape of the third molar tooth across all the age groups and not 
influenced by the grainy effect of the image. On the other hand, CV shows decreasing trend of segmentation 
accuracy from Group | (5-7 years old) to Group 6 (20-23 years old). In the first two groups (Group I and 
Group 2), the roots are not fully emerged, therefore CV is able to produce high segmentation accuracy. 
However, in other age groups, when the roots of the tooth start to grow, CV experienced difficulty to perform 
segmentation between the roots and gum, because their intensity contrast is low, which deteriorates the 
segmentation accuracy. Overall, average accuracy of DP-GICOV is 95.3%, and CV is 81.6%. 
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Figure 9. Results of classification accuracy of third molar teeth for persons aged from 5 to 23 years old 


4. CONCLUSION 

Digital imaging allows noninvasive forensic odontologist analyses to be performed. Out of the many 
teeth, third molar or wisdom tooth is the most reliable indicator for age estimation. In this paper, the third 
molar tooth is segmented using DP-GICOV and CV techniques. Samples acquired from different age groups 
are used as they show variability in terms of morphological structure as persons grow older. DP-GICOV 
technique produces better segmentation results than that of the CV technique. For the future works, focus will 
be given on the estimation of age that relate with the shape of the tooth that has been segmented. The 
outcomes of this work can potentially be used to identify an individual or a victim in dental forensic 
investigation. 
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